UCT-TP-289/12 



Chiral symmetry restoration and deconfinement in QCD at finite temperature 
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The light-quark correlator in the axial-vector channel is used, in conjunction with finite energy 
QCD sum rules at finite temperature, in order to (a) establish a relation between chiral-symmetry 
restoration and deconfinement, and (b) determine the temperature behavior of the ai(1260) width 
and coupling. Results indicate that deconfinement takes place at a slightly lower temperature than 
chiral-symmetry restoration, although this difference is not significant given the accuracy of the 
method. The behaviour of the ai(1260) parameters is consistent with quark-gluon deconfinement, 
as the width grows and the coupling decreases with increasing temperature. 
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I. INTRODUCTION 



The extension of the QCD sum rule method [l| to finite 
temperature was first proposed in [2j, and applied to 
the light-quark vector meson system, i.e. the p-meson 
channel. Additional theoretical support for the validity 
of such an extension was provided later in [3]. A quark- 
gluon deconfinement parameter was also introduced in 
[2| in the form of the squared energy threshold, sq(T), 
for the onset of perturbative QCD (PQCD) in hadronic 
spectral functions. Around this energy, and at zero tem- 
perature, the resonance peaks in the spectrum are cither 
no longer present or become very broad. The smooth 
hadronic spectral function thus approaches the PQCD 
regime. With increasing temperatures approaching the 
critical temperature for deconfinement, T c , one would 
expect hadrons to disappear from the spectral function 
which should then be described entirely by PQCD. The 
analysis of [2f indeed showed so(T) to be a monotonically 
decreasing function of the temperature, together with 
the coupling of the p-meson to the vector current. 
Since this analysis was performed in the zero-width 
approximation, an important dynamical feature also 
signaling deconfinement was overlooked, i.e. resonance 
broadening as proposed in [jM5|, and subsequently 
confirmed in other applications |fjj. 

A link between deconfinement and chiral-symmetry 
restoration using QCD sum rules in the axial-vector 
channel was first established in improved in Q, and 
recently updated and extended to finite density in Q. 
These analyses indicate that the temperature at which 
sq(T) vanishes is very close to that at which the quark 
condensate, or alternatively the pion decay constant 
f n (T) vanishes. Within the accuracy of the method 



these results imply that both phase transitions take 
place at roughly the same temperature. 
The analyses of 0-0] made use of the finite energy 
QCD sum rule (FESR) of the lowest dimension (d = 2) 
in the axial-vector channel, assuming full saturation of 
the hadronic spectral function by the pion pole. This 
assumption would not be entirely justified if one were to 
consider the subsequent two FESR of dimension d = 4 
and d = 6 in order to extract more information from the 
sum rules. In fact, already at T = one finds that the 
values of the condensates of dimension d = 4 and d = 6 
that follow from the second and third FESR are barely 
consistent with results obtained using experimental 
data [9j. This strongly suggests additional hadronic 
contributions, and in fact the data in this channel 
include not only the pion pole but also the ai(1260) 
resonance. A straightforward theoretical calculation 
confirms this to be the case. 



In this paper we reconsider the light-quark axial-vector 
channel using the first three FESR, together with an 
improved hadronic spectral function involving the pion 
pole as well as the ai(1260) resonance. This allows for a 
more realistic conclusion on the relation between chiral 
symmetry restoration and deconfinement. At the same 
time, it provides additional and valuable information 
on the temperature behaviour of the ai(1260) coupling 
and hadronic width. The results indicate that sq(T) 
vanishes at a critical temperature some 10% below that 
for chiral-symmetry restoration. Within the accuracy 
of the method this difference is not significant. The 
ai(1260) coupling initially increases with increasing 
T up to T/T c ~ 0.7, and then decreases sharply up 
to T c . The hadronic width of the cti(1260) remains 
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constant up to T/T c ~ 0.6, increasing sharply thereafter. 
This behaviour of the coupling and the width are fully 
consistent with a quark-gluon deconfinement scenario. 



II. FINITE ENERGY QCD SUM RULES AT 
T = 



We consider the correlator of light-quark axial- vector cur- 
rents 

n^(q 2 ) = i J dSe" 1 <0\T(A^(x) , 4(0))|0> 

= -g^U^q 2 ) +q fl q„n (q 2 ) , (1) 

where A^(x) =: d(x)j fl j5u(x) : is the (charged) axial- 
vector current, and = (u>,q) is the four- momentum 
carried by the current. The functions n 0) i(g 2 ) are free of 
kinematical singularities, an important property needed 
in writing dispersion relations and sum rules. Concen- 
trating on e.g. Tlo(q 2 ) and invoking the Operator Prod- 
uct Expansion (OPE) of current correlators at short dis- 
tances beyond perturbation theory, one of the two pillars 
of the QCD sum rule method, one has 

^IIo^Iqcd = C i+ £ C2N ^/ ] {6 2N {f)) , 

N=l ^ 

(2) 

where Q 2 = -q 2 , (6 2N (fi 2 )) = (0 |a 27V (a* 2 ) 1 0) , fi 2 is 
a renormalization scale, the Wilson coefficients Cn de- 
pend on the Lorentz indexes and quantum numbers of 
the currents, and on the local gauge invariant operators 
On built from the quark and gluon fields in the QCD 
Lagrangian. These operators are ordered by increasing 
dimensionality and the Wilson coefficients are calculable 
in PQCD. The unit operator above has dimension d = 
and CqI stands for the purely perturbative contribution. 
At T = the dimension d = 2 term in the OPE cannot 
be constructed from gauge invariant operators built from 
the quark and gluon fields of QCD. In addition, there is 
no evidence from such a term from analyses using the 
experimentally measured axial-vector spectral function 
Q. The dimension d — 4 term, a renormalization group 
invariant quantity, is given by 

Ci{6±) = ^(a s G 2 ) + 2w 2 (m u + m d )(qq), (3) 
6 

where the second term is negligible in comparison with 
the gluon condensate, and thus it will be ignored in the 
sequel. The leading power correction of dimension d = 6 
is the four-quark condensate, which in the vacuum satu- 
ration approximation [l| becomes 

C 6 (6 6 ) = 7 -^n 3 a s \(qq)\ 2 , (4) 



which has a very mild dependence on the renormaliza- 
tion scale. This approximation has no solid theoretical 
justification, other than its simplicity. Hence, there is no 
reliable way of estimating corrections, which in fact ap- 
pear to be rather large from comparisons between Eq. Q 
and direct determinations from data [9j. This poses no 
problem for the finite temperature analysis, where Eq.([4j 
is only used to normalize results at T = 0, and one is usu- 
ally interested in the behavior of ratios. 
The second pillar of the QCD sum rule technique is 
Cauchy's theorem in the complex squared energy s-plane, 
leading to the FESR (at leading order in PQCD) 

R (A, - 1} C 2N (0 2N ) = 4vr 2 / dss N - 1 -Imn ( S )| H AD 

Jo k 

S -l[l + 0(a s )} (A = l,2,-..). (5) 
The normalization of the correlator in PQCD is 
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FIG. 1: The ALEPH data in the axial-vector channel [Tlj . 
and in the resonance region together with the fit, Eq.© in 
the region relevant to the FESR. 

Imn ( S )| QCD = J- [1 + 0(a s (s))} . (6) 

4 7T 

At T = the radiative corrections above are known up to 
five-loop order, i.e. 0{a A s ), in PQCD. Higher dimensional 
condensates are poorly known [9( and thus will not be 
considered here. 

In the hadronic sector the spectral function involves the 
pion pole followed by the ai(1260) resonance 

Imn ( S )| H AD = 2^/2(5( s )+Imn ( S )| ai , (7) 

where U = 92.21 ± 0.14 MeV [To[ is the pion decay 
constant, the pion mass has been neglected, and a fit to 
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the ALEPH data [Tl| in the resonance region gives 



Im ITo(s)| 0l = C f ai exp 



Ml 



r 2 



(8) 



where M ai = 1.0891 GeV, T ai = 568.78 MeV, and 
C f ai — 0.048326. Using the first Weinberg sum rule as 
a rough estimate gives f ai = 0.073, and thus C — 0.662. 
Equation © is valid up to s = 1.2 GeV 2 , after which 
it becomes constant up to s ~ 1.5 GeV 2 . This fit to- 
gether with the ALEPH data is shown in Fig.l up to 
s ~ 1.5 GeV 2 (the FESR determine s = 1.44 GeV 2 ). 
The pion decay constant is related to the quark conden- 
sate through the Gell-Mann-Oakes-Renner relation 



2 / 2 M 2 = -{m u + m d )(0\uu + dd\0) 



(9) 



Chiral corrections to this relation are at the 5% level [12| , 
and at finite temperature deviations are negligible except 
very close to the critical temperature [l3j . 
The first three FESR can now be used in order to deter- 
mine the PQCD threshold sq, and the d = 4 and d — 6 
condensates. These results will subsequently be used to 
normalize all finite temperature results. The value of 
so obtained by saturating the hadronic spectral function 
with only the pion pole, and to leading order in PQCD, 
is s ~ 0.7 GeV 2 , as in [H, Q[H- This value increases 
substantially to a more realistic so = 1.15 GeV 2 once the 
ai(1260) contribution is taken into account, and becomes 
s = 1.44 GeV 2 with PQCD to five-loop order. Results 
for the leading condensates obtained from the FESR of 



dimension d = 4, 6 are in good agreement with direct 
determinations from data [9[ at the corresponding value 
of so, which is to be expected as Eq.([5]) provides a very 
good fit. 




FIG. 3: The gluon condensate C 4 (04}(TVC , 4(e>4}(O) as a 
function of T/T c from lattice QCD results [20]. Solid squares 
and circles correspond to 2 and 4 quark flavours, respectively, 
and error bars are the size of the points. The dotted curve is 
a fit to these data and the solid curve a smoothed fit. 
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FIG. 2: The quark condensate (qq)(T)/(qq)(0) = 
fl{T)/fl{0) as a function of TJT C in the chiral limit (m, = 
Mtt = 0) with T c = 197 MeV (solid curve), and for finite 
quark masses from a fit to lattice QCD results fl9l ] (dotted 
curve) . 



III. FINITE ENERGY QCD SUM RULES AT 
T / 



The extension of the QCD sum rule method to finite tem- 
perature implies a T-dependence in the OPE, as well as 
in the hadronic parameters. Regarding the former, both 
the Wilson coefficients as well as the vacuum condensates 
become temperature dependent. In particular, the strong 
coupling a s {Q 2 ,T) now depends on two scales, the ordi- 
nary QCD scale Aqcd associated with the momentum 
transfer, and the critical temperature scale T c associated 
with temperature. In principle, this poses no problems in 
the asymptotic freedom region and at very high tempera- 
tures, T >> T c , where PQCD can be applied. However, 
the QCD sum rule method approaches T c from below 
(starting from T — 0), so that the presence of this sec- 
ond scale is problematic. No satisfactory solution to this 
problem exists, so that analyses must be carried out at 
leading order in PQCD. This circumstance is of little con- 
sequence, since basically all hadronic parameters would 
hardly ever be measured with the same precision as at 
T = 0. At this order in PQCD there are two thermal 
corrections to Eq.®, namely one in the time- like region 
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(q 2 > 0), the so called annihilation term which in the 
static limit (q — > 0) is 



ImII+(u;,T) = — 

4 7T 



l-2n F 



(-) 

\2TJ 



(10) 



and one in the space-like region, the so-called scattering 
term, originating in a cut centered at the origin on the 
real axis in the complex energy u) = y^-planc of width 
— |q| < cj < |q| [2]. In the static limit this is given by 



Imn (w,T) = ^ S(uj 2 ) J yn F (J^ dy 



T 2 S(co 2 ) 



(11) 



where np(z) = l/(l + e z ) is the Fermi thermal function, 
and the chiral limit was assumed. 
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FIG. 4: Fhe continuum threshold so(T)/so(0) signaling de- 
confinement (solid curve) as a function of T/T c , together with 
fl{T)/fl{0) = (qq){T)/(qq)(0) signaling chiral-symmetry 
restoration (dotted curve). 



couplings, and widths become T-dependent. Hadroni- 
cally stable particles, e.g. the pion, with T(0) = develop 
a width , although this effect is far less pronounced than 
in the cases where T(0) ^ 0. The important parameters 
signaling deconfinement are the hadronic width and cou- 
pling, but not the mass. In fact, the latter is just the real 
part of the hadron propagator in the complex squared en- 
ergy plane, while the width is its imaginary part. A van- 
ishing mass at T — T c would not signal deconfinement, 
unless the width diverges at such a temperature. But 
then the value of the mass becomes irrelevant, it could 
just as well retain its zero temperature value, or increase. 
This is actually what QCD sum rule analyses show [f|, 
i.e. the hadronic mass is essentially constant in a wide 
range of temperatures, increasing or decreasing slightly 
very close to T c , depending on the channel, and the width 
diverges at T = T c . A notable exception are the scalar, 
pseudoscalar, and vector charm-anti-charm states which 
survive beyond T c 17]. Finally, there is in principle a 
hadronic counterpart to Eq. fTTj) originating in current- 
pion scattering. However, in the axial-vector channel this 
term is loop suppressed, as the current only couples to 
an odd number of pions. 

The temperature behaviour of the quark conden- 
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At finite temperature there are in principle additional 
contributions to the OPE, Eq.([2]), in the form of nondi- 
agonal (Lorentz noninvariant) condensates. In the case of 
nongluonic operators these are highly suppressed [|| , 
so that they can be safely ignored. A gluonic twist-two 
term in the OPE was considered in [lj| , but it is at least 
two orders of magnitude smaller than the standard gluon 
condensate at the temperatures considered here. At di- 
mension d = 2 there is evidence for the presence of a non 
gauge invariant condensate at high temperatures [lj|. 
However, at the temperatures explored in the present 
analysis this condensate can be safely neglected. 
In the hadronic sector and at finite temperature, masses, 



FIG. 5: The hadronic width of the ai(1260) resonance 
r ai (T)/r ai (0) as a function of T/T c . 

sate, equivalently f 2 , is shown in Fig. 2 in the chiral 
limit (solid curve) as determined in the framework of 
the Schwinger-Dyson equation (l8| . and for finite quark 
masses from a fit to lattice QCD results (dotted curve) 
0, (HI. The critical temperature is T c = 197 MeV. In 
the sequel we concentrate on the chiral limit as we find 
that the FESR have solutions only up to T ~ 0.9 T c , a 
region where the quark condensate is essentially unique. 
The gluon condensate is shown in Fig. 3 from a fit to 
lattice QCD determinations (dotted curve) [2(| , together 
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with a smoothed fit, both adjusted to T c = 197 MeV. 
This smooth fit is needed to avoid instability in the FESR 
due to the sharp break in the lattice QCD curve. 
Making use of the above information the first three ther- 
mal FESR become 



^fl{T) = |vr 2 T 2 + J da 



2n F — 



4tH 



«o(T) 



1 



ds - lmU (s,T)\ ai 



C 4 (0 4 )(T) = 4vr 2 / dss-Imn (s)| ai 

Jo T 



ds . 



1 - 2n F — 



(12) 



,(13) 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

T/Tc 



fMT) i 

C 6 (0 6 )(T) = 4tt 2 / dss 2 -Imn (s)| Ql 

JO 7T 



*o(r) 



(is s 



2n F — 



(14) 



These equations determine the continuum threshold 
so(T), the coupling of the ai(1260) to the axial- vector 
current, f ai (T), and its width T ai (T), using as input 
the thermal quark condensate (or / 2 (T)), the thermal 
d = 4, 6 condensates, and assuming the ai(1260) mass 
to be temperature independent, as evidenced by results 
m many different channels @, [H|. 



IV. RESULTS AND CONCLUSIONS 



FIG. 6: The coupling of the oi (1260) resonance f ai (T) // 01 (0) 
as a function of T/T c . 



method. The behaviour of the width is shown in FigfSJ 
and that of the coupling in Figf^J The rise of the 
width, and the fall of the coupling are indicative of a 
transition to a quark- gluon deconfmed phase at T = T c , 
and provide additional support for the interpretation 
of sq(T) as a phenomenological order parameter for 
quark-gluon deconhnement. It should be stressed that 
resonance broadening as obtained here is the result of 
an interplay between QCD and hadronic information. 
Hence, it is directly related to quark-gluon deconhne- 
ment, in contrast to width results at finite T from purely 
hadronic models such as e.g. the sigma model [22], which 
reflect a purely hadronic (absorption) effect in a medium. 



The FESR have solutions for the three parameters, 
s a (T), f ai (T), and r oi (T), up to T ~ (0.85 - 0.90) T c , 
a temperature at which sq(T) reaches its minimum. An 
inspection of Figf2] shows that at such temperatures 
the thermal quark condensate, or equivalently f n (T), is 
independent of whether the chiral limit (massless quarks) 
is assumed or not. A short extrapolation to T — T c is 
to be understood for all results in the sequel. An in- 
spection of Eq. (jT2l shows that disregarding the ai(1260) 
contribution, so(T) would vanish at a lower critical 
temperature than f n (T) (or (qq)(T)). In fact, making 
the very rough approximation of neglecting the thermal 
factor iif(t/s/2T) in the second term on the r.h.s. of 
Eq.(HU leads to s (T) ~ 8tt 2 / 2 (T) - (4/3)tt 2 T 2 . This 
feature remains valid even after including the ai(1260) 
in the FESR, as shown in FigfJJ corresponding to the 
solution for s (T) using all three FESR. . In any case, 
this 10% difference is well within the accuracy of the 
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